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We propose a mechanism for testing the theory of collapse models such as continuous spontaneous 
localization (CSL) by examining the parametric heating rate of a trapped nanosphere. The random 
localizations of the center-of-mass for a given particle predicted by the CSL model can be understood 
as a stochastic force embodying a source of heating for the nanosphere. We show that by utilising a 
Paul trap to levitate the particle and optical cooling, it is possible to reduce environmental decoher¬ 
ence to such a level that CSL dominates the dynamics and contributes the main source of heating. We 
show that this approach allows measurements to be made on the timescale of seconds, and that the 
free parameter Acsi which characterises the model ought to be testable to values as low as 10“^^ Hz. 


Dynamical reduction models - better known as col¬ 
lapse theories - seek to resolve the measurement prob¬ 
lem by inserting a non-linear stochastic term in the 
Schrodinger equation (SE). This would account for gen¬ 
uine collapses of superposition states. In these theories, 
'localisation' events occur at a frequency scaling with 
the mass of the system at hand. These are frmdamen- 
tally different to environmental decoherence (l|], and are 
invoked as the origin of wavefunction collapse. Such 
modifications to the SE aim to provide a theory capable 
of describing phenomena at all scales, and are designed 
to reproduce conventional quantum mechanics (CQM) 
when dealing with small masses, and classical mechan¬ 
ics at the macroscopic scale. 

One of the most celebrated models of dynamical re¬ 
duction is the Continuous Spontaneous Localization 
(CSL) model i]. It is characterised by two parameters: 
a length Tc, and a frequency Acsi- The former provides 
a length scale above which reduction effects would be 
relevant, the latter embodies the rate at which spatial 
superpositions of a single nucleon separated by a dis¬ 
tance greater than Vc would collapse. While Vc is gener¬ 
ally taken to be ~ 100 nm, the value of Acsi is the subject 
of uncertainties lii] and is currently taken to span a 
range from 10“^® Hz 10] to 10“® Hz |@]. A value of 
10-8±2 ]3ggj-( proposed @], based on the process 

of image formation on photographic film. The heating 
rate of ultracold atoms was used to set a value of 10“^ 
Hz 0, while the maximum allowable heating rate of 
the intergalactic medium seems to be compatible with 
10-io±2 ji]. 

The effects of localization are mathematically very 
similar to those of decoherence |0], meaning that any 
experiment built to search for a signature of such a col¬ 
lapse mechanism must minimize the effects of decoher¬ 
ence as much as possible so as to better distinguish the 
hallmark of the former from that of the latter. However, 
for a given object, both its rate of localization and de¬ 
coherence will increase proportionally to mass. To test 
collapse we must study objects large enough to have an 


appreciable localization rate, yet small enough that de¬ 
coherence does not dominate the d5mamics. The scale at 
which this becomes possible is the so-called mesoscopic 
one, the liminal scale between the well established quan¬ 
tum and classical regimes. Recently, Nimmrichter et al. 
have shown iB that beyond a certain size collapse ef¬ 
fects have a sub-linear scaling with size, a result which 
we corroborate here, and which clearly identifies the 
regime of interest. 

Lately there have been a considerable number of pro¬ 
posals to test collapse theories, and the challenge of re¬ 
ducing decoherence manifests itself differently depend¬ 
ing on the scheme iiMi. Unfortunately many pro¬ 
posals do not include detailed analyses of conventional 
environmental noise, and as such, it is not clear what 
range of values of Acsi they could probe. Matter-wave 
interferometric methods are an attractive means for such 
tests with a well established range of testable values of 
-^csl B. Such settings, however, require the preparation 
and detection of many identical particles, which makes 
implementation challenging. Optomechanical propos¬ 
als 10] are attractive in comparison, because they in prin¬ 
ciple require only a line shape measurement of the light 
scattered by an intra-cavity macroscopic oscillator, and 
do not require ground state cooling. However, as such 
line shape would be narrow (order of /iHz) very long 
times (in excess of months) would be required for the 
measurement. A recent proposal has been made based 
on dynamical decoupling 12011 . which is promising but 
may also be constrained by long testing times. 

The concept of utilising the energy gain of a harmonic 
oscillator to test CSL was first suggeste d by Adler ll2l|] 
and more recently re-examined in Refs. i22f. 12311 . Whilst 
this method is relatively straightforward and therefore 
attractive for implementation, the ability to test for CSL 
is very dependent on a detailed and realistic inclusion of 
conventional decoherence mechanisms. In order for the 
classical approach put forward, for instance, in Ref. 11231] 
to be effective, the collapse mechanism must be able to 
induce the excitation of a rather substantial number of 
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FIG. 1. Schematic of the experiment, in which the particle is 
levitated by the electric field of the Paul trap, and cooled by 
the optical cavity. 


thermal phonons in a given oscillator. This is, in general, 
not the case for a large range of values of the parameters 
that characterize a given collapse model and such ap¬ 
proaches can only assess the largest of their conjectured 
values. 

In this paper, we propose an experimentally viable 
way to explore CSL on the mesoscale by utilising a 
cavity-cooled, single-charged nanosphere trapped in a 
Paul trap j^. Measurements can be made with a sin¬ 
gle trapped particle in less than 100 seconds and, un¬ 
der optimal conditions, we find fhis scheme capable of 
probing Acsi to values as low as 10“^^ Hz, thus going 
significantly beyond the literature reported so far. Most 
importantly, our protocol allows for the discrimination 
between collapse effects and mis-characterised conven¬ 
tional noise ones- a challenge which has to our knowl¬ 
edge remained unaddressed so far. 

We explore the possibility for an optomechanical test 
of a form of CSL described in Ref. 0, ll^ . in which 
the effects of spontaneous localization are modelled as 
a delta-correlated stochastic noise source wt- This ap¬ 
proach is valid when the scale of spatial-superposition 
separations is less than Tc- The noise term wt will oc¬ 
cur in the dynamical equations of the system as an extra 
Langevin force EM- Its effect on the dynamics of a me¬ 
chanical oscillator would depend on the size and den¬ 
sity of the object collapsing, and the two parameters Vc 
and Acsi characterizing the model. Conveniently, we can 
represent the effects of the localization process via a dif¬ 
fusion operator characterised by the coefficient Dcsi and 
appearing in the master equation describing the particle 
in the same way as conventional heating sources. The 
diffusion coefficienf takes the form 


Hcsl 


Ti Acsi 
mwm rj. 


( 1 ) 


where a is a geometry-dependent factor, which for a 
sphere is given by 


m 

mo 
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Here mo = 1 amu and R is the sphere's radius. 

The Protocol- We now describe the scheme that we pro¬ 
pose to test the CSL model. We levitate a charged 
nanosphere in a hybrid trap consisting of a Paul trap 
and an optical cavity, and use them to cool its motion 
to a temperature corresponding to a low occupation 
number. We then turn off the optical field (and hence 
the cooling) and let the dynamics evolve for a certain 
amount of time before measuring the energy of the os¬ 
cillator again. A model including the effects of CSL pre¬ 
dicts it will have heated more than one would expect 
due to conventional noise sources alone. If the mea¬ 
sured energy matches that predicted by conventional 
noise sources, we will have provided evidence against 
CSL to within a certain range of Acsi, whereas a higher 
measured energy would indicate some other dynamics 
at play, in favour of collapse theories. 

We divide the procedure into two phases: a cooling 
phase, and free evolufion. For fhe purpose of testing 
CSL, it is the second phase that is important. In this 
period of free evolution the nanosphere is levitated us¬ 
ing a single electric potential, which could be provided 
by a number of generic trap architectures. The mecha¬ 
nism of cavity cooling for nanoscale objects is well es¬ 
tablished 10-10 Hi 12 ^^. 1^ . and relies on having the 
particle sit in two potential wells: one of which traps, 
and one (or both) of which cool. Though these poten¬ 
tials are traditionally provided by an optical cavity pop¬ 
ulated with two distinct optical fields. Ref. Hi shows 
that cooling is possible using a Paul trap in conjunction 
with a single-mode optical cavity. While both these sys¬ 
tems are required for the cooling phase, the Paul trap 
alone suffices to levitate it. 

This is appealing, because in the low pressure sce¬ 
nario of a particle levitated solely by optical fields the 
dominant source of heating is the scattering of cavity 
photons B- By using a hybrid trap we can cool the 
particle to a desired temperature and then turn off the 
optical field completely, leaving the particle suspended 
in the Paul trap alone. This ability to turn off fhe opti¬ 
cal field without 'losing' the particle means that we can 
do away with what would otherwise be the dominant 
cause of heating - optical scattering - and thus greatly 
reduce the conventional heating sources that would oth¬ 
erwise mask the CSL effects. 

We emphasize that although we require cooling, we 
do not need to achieve the ground-state energy. Indeed, 
the simple comparison between the initial phonon num¬ 
ber no and the final one nf affer fhe period of free evo¬ 
lufion will give us information on the heating rate B- 

The period of free evolution is governed by the 
Hamiltonian H = Hq + H', where Ho = a, aJm is 

the secular frequency of the Paul trap, , a are the cre¬ 
ation and annihilation operators for fhe centre of mass 
motion of the sphere respectively, and H' represents the 
interaction between system and environment. We can 
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then solve the master equation p = —{i/h)[H,p] for the 
oscillator 10. The forms of coupling to noise sources in 
H' determine their effects on the master equation 0. 
We have explored each noise source in detail, examin¬ 
ing collisions with the background gas, blackbody ra¬ 
diation, acoustic noise affecting the trap, Johnson and 
patch potential noise from the electrodes, micromotion 
from the trap's driving frequency, and anisotropy of the 
sphere 0- We group these noise sources as momen¬ 
tum diffusion, occurring at rate Udiff/ position diffusion 
at rate Upos/ and momentum dissipation at rate T, thus 
getting the dynamical equation 14011 

2 

p = - J[Wo,p] - E [X,,p\] - r[Q., {A, p}] (2) 

f=l 

where X = {Qz,Pz) is the vector of quadratures of 
the nanosphere Qz = a + and Pz = i{a^ — a), and 
D = (Z?diff, Tlpos). Our analysis of the various noise 
sources, including the possible effects of anisofropy of 
fhe nanosphere (and the consequent non-uniform dis- 
fribution of the charge) finds all but gas collision and 
blackbody radiation to be negligible llTI] , giving us T = 

(Tgas Tbb,e “t” Tbb,a)/d aud -f^diff ~ Pgas T Phh T Pcs\/ 

where D^si is given in Eq. 1(1) and 

^g^sPenv ^ ^B {'Jhh,ePint qbb,aTenv) 

= 2nuj^ ’ ■ 

(3) 

Here, 7gas, 7bb,e, 7bb,a are the damping constants related 
to gas collisions, blackbody emission and blackbody ab¬ 
sorption respectively. The environmental temperature is 
Tenv, ujm and m is the mass of the nanosphere. 

Heating Rate of a Trapped Particle- The mean occupation 
number (phonon number) of the trapped particle as a 
function of time (n)t is determined via Eq. (|2). Assum¬ 
ing a thermal state ESS, the expression for the 
rate of change of the average phonon number simplifies 
fo {n{t)) = —T{n(t)) -t- Ddis, which has the solution 

(„(i)) = (n. (4) 

where ng is the initial average number. In Eig. |2]we plot 
the expected mean phonon number over the first sec¬ 
ond of free evolution, starting from an initial number 
of no = 50, showing the heating when the CSL mech¬ 
anism is included (solid blue line), and when it is not 
(dashed orange line). In our case, the diffusive terms in 
Eq. (|U dominate over the dissipative, resulting in an ex¬ 
pression that is approximately linear, as seen in Eig. |2l 
The inclusion of the CSL mechanism results in a heating 
rate of the nanosphere motion of about 2500 phonons/s, 
which is in stark contrast with the Acsi = 0 case, where 
only ^ 350 phonons/s are achieved for the parameters 
used in our simulations. 


Differentiating CSL from Conventional Noise- If a final oc¬ 
cupation number is recorded that agrees with the model 
of Acsi = 0, we interpret it as falsifying a certain range of 
Acsb However, if we measure a higher phonon number, 
we can infer that some extra heating process is present, 
possibly collapse effects. However, an objection could 
be made that the increased heating would more likely 
result from mis-characterisrng the environmental noise 
sources present. 

Therefore, an important requirement is how to cor¬ 
rectly identify CSL, and distinguish it from other noise 
sources. Indeed, this problem is generic to any test of 
collapse theories. We can address this problem by mon¬ 
itoring the conventional noise sources, such as black¬ 
body radiation and gas collisions, by varying the associ¬ 
ated parameters and determining the effect they have on 
the heating rate. For example, if the heating is strongly 
dominated by CSL, then varying the pressure will have 
little effect over some range, while a gas-collision dom¬ 
inated process results in an almost linear relation be¬ 
tween pressure and heating rate. 

Fig. |3] summarises such behaviours, and displays the 
trend followed by the mean phonon number after one 
second when varying different parameters for the cases 
of Acsi = 10“® Hz and for Acsi = 0. In Fig. |3] (a), we 
examine the effect of varying the background pressure 
when we include or exclude CSL. The difference in re¬ 
sponse between the theories is instructive. Without CSL 
the effect of an increasing pressure can be seen across the 
whole range, whereas for CSL there is a region of im¬ 
munity where CSL dominates the dynamics. Likewise 
for the internal temperature depicted in Fig. |3] (b): the 
point at which this significantly influences the heating 
is different for the two theories. The response to a vary¬ 
ing mechanical frequency also takes a different shape, 
as seen in Fig. |3] (c). Most interesting is the effect of a 



FIG. 2. (Color online) Expected phonon number {n{t)) over 
one second with and without collapse effects. We have used 
no = 50, R — too nm, a pressure ofl0“^^ mbar, ujrn = 5 kHz, 
p = 2300 kg/m®, Tenv = 4 K, and an internal temperature of 
65 K. The solid (dashed) line is for a CSL mechanism charac¬ 
terized by Acsi = 10“®, (0) Hz, i.e. the Adler value (no CSL 
mechanism). 
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varying radius shown in Fig. |3](d). Our findings agree 
with those of Ref. M- objects must be large enough to 
have an appreciable collapse rate, but small enough that 
decoherence does not dominate the d5mamics in order 
for us to observe collapse effects. Such request is met for 
R G [10, lOOjnm, roughly. For all of these we have com¬ 
pared the Adler value for Acsi = 10“® Hz with the case of 
no collapse effects. However the same differences per¬ 
sist for any chosen non-zero value of Acsi, though being 
more pronounced for higher collapse rates. 

In order to find the best conditions for testing CSL 
it is necessary to numerically optimise all the param¬ 
eters simultaneously. In Fig. IH we show the range of 
Acsi which can be probed for an illustrative set of exper¬ 
imental conditions. We see that, as one would expect, 
the testable range depends ultimately upon the condi¬ 
tions that can be achieved, most relevantly the mini¬ 
mum value of both environmental pressure and inter¬ 
nal temperature. The value of the initial phonon num¬ 
ber and evolution time also play a significant role: af¬ 
ter enough time, the phonon ratio will tend to the heat¬ 
ing rate ratio as (n(oo))csi/(n(cx)))cqm -t (h)csi/(h)cqm, 
and the time required to approximate this depends on 
the initial phonon number and heating rates. A lower 
initial phonon number or longer evolution times would 
promote each plot in Fig. IH as a given set of parameters 
would be capable of probing a lower value of Acsi- 
Experimental Feasibility. - Charged silica particles of 200 
nm have been trapped and cooled to milliKelvin tem¬ 
peratures in a Paul trap using cavity cooling to mil¬ 
liKelvin temperatures jll] thus demonstrating the key 
experimental components required for our proposal. 
Pressures down to the 10“^^ mbar range and internal 

(a) (b) 






FIG. 3. (Color online) Phonon number after Is as predicted 
by with (blue line) and without (dashed orange line) collapse 
effects as we vary the pressure in panel (a), bulk temperature 
in panel (b), mechanical frequency in panel (c), and radius of 
the sphere in panel (d). Except where a parameter is under 
investigation, we have used the same values as in Fig.lD 



FIG. 4. (Color online) Testable range of Acsi as a function 
of pressure, for a set of internal temperatures. The pressure 
shown is where (n(100s))csi/(n(100s))cqm > 1-2 after 100 s 
with initial phonon number of no = 50. The sphere radius 
and mechanical frequency have been optimised for each indi¬ 
vidual data point's temperature and pressure. 

particle temperatures in the 10 K range can be obtained 
for this setup using standard cryopumps. An important 
component of these experiments is the measurement of 
the oscillator energy following parametric heating by 
CSL and conventional noise. These can be carried out 
by performing a homodyne measurement of position as 
a function of time which allows the determination of 
mean occupation number. Importantly, the light-cavity 
detuning can be non-perturbative during the measure¬ 
ment such that it neither cools nor heats the particle 
during this time. This ability to control the cooling rate 
also allows us to avoid giving strong kicks to the parti¬ 
cle when we turn off the optical field after the particle is 
cooled to the desired energy. 

Concluding Remarks.- We have shown that the paramet¬ 
ric heating rate of a trapped nanosphere provides a 
viable mean of testing CSL. Central to the success of 
this scheme is the minimization of all sources of en¬ 
vironmental decoherence. The two-stage 'cool and re¬ 
lease' protocol that we have illustrated allows us to ex¬ 
ploit optical cooling while avoiding problematic scatter¬ 
ing noise that would otherwise dominate the dynam¬ 
ics. We remark that, owing to established techniques 
and common experimental settings, the experiment can 
be performed rapidly, and is repeatable upon a given 
nanosphere. Remarkably, the central e^erimental set¬ 
up has already been demonstrated Our scheme 

can be used to distinguish CSL from other noise sources 
- an essential condition for inferring the existence of 
collapse effects from an experiment. Further, we have 
shown that the parameter range of CSL that is testable 
using our proposal is broad, and can readily be expected 
to probe Acsi = 10^° (Acsi = 10“^^ ) using a background 
pressure of 10“^^ mbar mbar) and an internal 

temperature of 60 K (20 K). Based on state of the art, 
values as low as Acsi = 10“® could be tested imminently, 
and Acsi = 10“^^ plausibly in the next few years. 
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COOLING THE SPHERE 


Ref. is^ give a comprehensive mechanism for cooling the motion of a sphere in a hybrid type trap. Though 
they achieve an experimental temperature of 10 K, improved to millikelvin in llsil] . the same technique should in 
principle be able to achieve ground state cooling and is limited primarily by the firmesse of the mirrors. When 
translating their results to our scenario however, we must be a little careful. Along the axis of interest, there are three 
motional frequencies fo be considered in the hybrid scheme: utd, which is the frequency of the AC voltage applied to 
the electrodes to keep the trap stable; Wg (or secular frequency), which describes the oscillation of the sphere in the 
electric potential; Wc, the frequency of fhe sphere's motion inside a well of the cavity potential. For our period of free 
evolution, we are considering the motion inside the electric potential Wg = ujm, whereas th e goal of Ref. was to 
cool ujc- The steady state phonon number achievable through cooling is given by IS12LIS13ll 

»r _ f Ksc \ Fsc Tothers /q.. .. 


where we have used N to refer to the phonon number during the cooling period and n for the phonon number dur¬ 
ing free evolution. Tothers refers to heating processes other than optical scattering, which are all negligible compare 
to rsr, and Ksc is the rate of photon loss from the cavity due to scattering, and the time averaged cooling rate is given 

byQ 


T_=g^k 


_(A-CCc)2 + f (A+CCc)2 + f_ 
in which k is the wavenumber of the laser, and the optomechanical coupling is given by 


(S2) 


9^ = 


1-Jo(4fcAd).,2 2 /3Fe.-l 

■rik a„ 


2mu}c 


214 e + 2 


UJl 


(S3) 


where is the number of photons, of frequency oji circulating in the cavity, whose mode volume is 14- Here, er is 
the relative dielectric permittivity, Jo is the zeroth Bessel function of the first kind and is the amplitude of the AC 
drive voltage. For a firmesse of 10®, grormd state cooling within the optical potential should be achievable using an 
input laser power of 10“^ W. Using {N)(t) = (iV)therme“^“* + N^s (l — e'"-*) where (iV)therm is the mean phonon 
number the oscillator would have if thermalised with its environment, we can see that {N){t) « {N)ss within about 
a microsecond, so there is no time impediment to achieving the cold states we require. 

Cooling to a given phonon number in the optical potential does not translate directly into that phonon number 
being occupied in the Paul trap once the cavity field is turned off, due to the different trap frequencies. Supposing 
fhe antinodes of the optical well and Paul trap are perfectly aligned, when we turn the optical potential off the sphere 
should have uq = NssUic/ujs, where uq is the phonon number in the Paul trap immediately after the cavity is switched 
off and and Nss is the phonon number in the optical trap immediately before. However, if we take into account some 
displacement between the two trap centers 6x, then we would expect the that no > NssUJc/ujs + mujsSU{2K), where 
the inequality accounts for the fact that it is also possible for the optical field to impart momentum to the sphere as 
it is turned off. The ability for this displacement to impart phonons requires that we ensure ^ 0.5nm. 

When we consider these factors, it becomes clear that our predictive knowledge of fhe phonon number at the 
begirming of the free evolution no is far from perfect. This is important of course, because it is by comparing no with 
nf that we hope to learn anything. However, it is possible for us to gain experimental, not just predictive knowledge 
of no- A 'dry run' is possible, in which rather than turning the cavity field completely off we reduce it to such a low 
power that it no longer traps the particle, but d oes re main coupled to the particle's position. As such we can use it 
to measure no via the techniques laid out in Ref. iS2ll] . This information can be used heuristically to better align trap 
centres, and also to build up a statistical picture of no, providing a benchmark for the actual experiment in which the 
cavity field is turned off complefely. 
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THE CSL MECHANISM AS A SOURCE OE NOISE 


The premise of the experiment is that since collapse effects can be thought of as physically, if not ontologically, 
equivalent to those of decoherenc^il a search for collapse is in some ways equivalent to a search for a new source 
of decoherence. More specifically, in our case we are searching for a Brownian noise source which heats a harmonic 
oscillator . He re we demonstrate that frequent random localisations can be modelled as Brownian noise, drawing on 
Refs. |S8, SI9II . The modified master equation is given by p = (U + Ccs\)p, where 


^cs\P - 


'^csl 


TrS/V^simg 


(TX 


M(x)pM(x) - -{p,m2(x)} 


(S4) 


and M(x) = exp[- 




:pr- 

' csl 


Here r„ = r + + Ar„ stands for the relative coordinate of particle n in relation 


to the centre of mass (C.O.M) coordinate r, with r° and Ar„ describing the motion of the particles bound within the 
solid. M (x) is a Gaussian averaged mass density function of the N particle system. As shown in Sec. 8.2 of Ref. js^l, 
by taking the range of motion for each constituent particle to be much smaller than Tc, it is possible to effectively 
separate the internal dynamics from those affecting the centre-of-mass (COM), meaning that we can re-write the CSL 
master equation in terms of the C.O.M alone, neglecting the internal dynamics. We thus have 

AcsiP = - p), (S5) 


where p(k) is the Fourier transform of the object's mass density. If we expand the exponentials in lISSll to first 
order an d the center of mass oscillations are taken to be much smaller than r^, we can simplify the CSL diffusion 
operator IS19ll to 


Acsip = -Dcsl[x, [x,p\]/h'^ = -oAcsl 


i,p]]/n‘^ 


(S6) 


where a is a geometric factor describing the shape and density of the object, which for a sphere (following again 
Ref. Q) is given by 


) :< 




■ 


(S7) 


We can now incorporate Dcsi into the master equation describing a macroscopic oscillator in the same way as any 
other diffusive noise source. 


Dissipation and Diffusion in the CSL Mechanism 

Most noise sources will cause momentum diffusion and momentum dissipation. Typical Brownian noise such as 
background gas collisions and blackbody radiation take this form, resulting in energy absorption (at rate D^ig) from 
the environment surrounding the nanosphere, as well as deposition of excitations into it (contributing to momentum 
damping at rate T). Other noise sources consist of one way mechanisms that result in unidirectional feeding of energy 
to the nanosphere (thus contributing to the value of Hdis) and thus heat it up without dampin g its motion. A relevant 
source of this form of noise would be the scattering of photons from the cavity field mode iS13ll : the oscillator can 
scatter these photons out of the cavity and gain energy without giving up anything to create them. 

CSL, as a developing theory, has been confronted with various problems, the question of energy conservation 
amongst them. As we've stated, the mechanism of localization in CSL arises via a coupling to a noise field. Through 
this constant one-way coup ling energy is transferred, and ultimately violates energy conservation. This issue was 
recently addressed in IS20ll . in which a new modification to the Schrodinger equation is introduced, resulting in 
a dissipative effect which remedies this problem and restores energy conservation. The rate of dissipation for a 
nanosphere of i? = lOOnm is found to be ^ 10“^ when using the Adler value for Acsu which is negligible for 
our considerations, and hence we leave it out of our analysis. As such, we treat noise from collapse as a solely 
contributing to momentum diffusion, and not to dissipation - making it analogous to optical scattering. 


^By this we mean that the effects of localisations are, at the level of the density matrix, indistinguishable from those of decoherence. Collapse willgppg^j. jg g] 
destroy coherences and localise a wavefunction. However, the meaning of this is profoundly different to usual: the system is being driven into a 
series of definitive states according to certain probabilities, it is not just being reduced to superposition states who have lost their coherence, and 
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CONVENTIONAL NOISE SOURCES 


Collisions With the Background Gas 


Collisions with background gas are typically the dominant source of noise for trapped nanospheres (following 
optical scattering). Since the nanosphere can both give and receive momentum to and from the background gas, it 
contributes both to Z?diff and F. The contribution to momentum diffusion is given by 


_ Tenv 

flUJm 


(S8) 


where Tenv is the temperature of the environment and 7 g is given in Ref. ISlOl] as jg = ■ Here, Vg is fhe mean 

velocity of the gas particles, d is the density of the sphere and P is the background gas pressure. The contribution to 
momentum dissipation is simply Fgas = 7g/4. 


Blackhody Radiation 


Blackbody radiation also contributes to both diffusion and dissipation, being a co uplin g to a bath with infinite 
modes at the temperature of the environment. The damping coefficients given in Ref. ISIDII are 


7bb,z 


27r4 

63 c^hdujrn e + 2 


(S9) 


where i represents either emissive or absorptive processes. We have Ti = Tenv for absorptive radiation and T^ = Tint 
for emissive. This gives us the diffusive term Dhh = ( 7 bb,eTint + 7bb,aTenv)fcs/(27iwm) and the dissipative one Fbb = 

(ybb.e ybbTa)/"!. 

We find thaf emissive radiation has a far sfronger effecf than absorptive. We can characterise the internal temper¬ 
ature which is primarily det ermi ned by the intensity of the cooling laser. The expression for the steady-state bulk 
temperature is given in Ref. ISllh as 


Tbulk = ( 




362 


euAlaser (ci "I" 2)^ -|- i 


1/4 


(SIO) 


where Iq is the intensity of the cooling laser at its waist and Aiaser is its wavelength, e is the emissivity of sphere, 
e = ei +162 is its complex dielectric constant, cr is the Stefan-Boltzman constant. 


Electric Eield Noise 

There are three mechanisms which will contribute to noise in the electric field experienced by the nanoparticle: 
Johnson noise , caused by thermal fluctuations in the electrodes; patch potential noise, thought to be caused by 
adatoms lISlSll on the surface of the electrodes; and acoustic vibrations in the lab, which will shake the electrodes 
relative to the trapped particle. Though the origins of these forms of noise are fairly independent, and their spectra 
will be different, all three will contribute to the heating of the particle and in practice it will be diffi cult to distinguish 
between them. They can all be modelled as a variation in the eff ective sp ring constant of the trap lS14ll . In ion traps, 
electric field noise is often dominant over other forms of heating ISlfillSlTjl . The heating rate due to patch fluctuations 
is given by 


patch 


Amhuj 




(Sll) 


where Se^ (w) = A(r)i?(w) is the power spectrum of the electric field noise, determined a geometry factor A(r) and 
the noise spectrum R^uj). A discussion of A(r) can be foimd in ISlSll . Whilst the details vary significantly between 
electrode archit ectur es, a trend is shared in that A(f^ oc d~°‘, where a ^ i. We can describe the noise spectrum of 
the adatoms as ISlSll R{uj) = aSfj_{uj), where a gives the density of adato ms on a surface and Sg, is the spectrum 
of the fluctuating dipole of a single adatom. By utilising the framework of lISlSll to describe the geometry of a trap 
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architecture, and that of BSlSil to describe the patch potentials which form on the electrode surfaces, a detailed picture 
can be built of the noise one would expect from any reasonable shape of trap, made of any metal. However, since we 
are not proposing a specific electrode architecture, it makes more sense to translate experimentally verified heating 
rates from other experiments into our own setting. Ref. IS16l] finds heating rates ~ 10 Hz for a single ion with 

a mechanical frequency ui'^ in the range of MHz. Assuming a similar trap geometry, we can see via Eq. JSllb that 

we would expect a heating rate of HEfieW = -PEfieid ~ 10“^, which is negligible compared to th e oth er noise 

sources that have been identified. Note that this prediction is based on experimental data gathered in llS16l] . and as 
such must reflect the presence of all three forms of electric field noise. 


Anisotropy of the Sphere 

The question of sphere anisotropy is a recognized problem when considering the dynamics of a trapped 
nanosphere, and is addressed in the supplementary information of lISlOll . In their analysis, the anisotropy presents 
a problem because it will cause an increased polarizability. This problem is resolved in our case, since there will be 
no optical field over the period of interest. The potential problem arising from anisotropy for us will be the corre¬ 
sponding anisotropy of the charge distributed over the surface, though we remark that this distribution may also be 
anisotropic even for an ideal sphere. This question of charge distribution is addressed below. 


Rotational Dynamics 

The anisotropy of the charge distribution on the sphere will cause a dipole moment, and this dipole moment will 
experience a torque due to the electric fields, and further it will affect the COM motion. The equations of motion for 
the COM are given by 


mr — gE(r, t) — V[/x • E(r, t)] = 0 (S12) 

where r = {x, y, z) is the position of the nanosphere in three dimensions, /r is its dipole moment, and E(r,t) is the 
electric field at point r at time t. Clearly we can see a coupling between the orientation of the dipole moment and 
the COM motion. This could potentially constitute a heating source. In order to understand what this heating rate 
amounts to, we need to model the rotational dynamics of the dipole together with the translational motion of the 
sphere, and see what the effect of the dipole moment is upon the COM motion. 

This situation is very similar to that of diatomic molecules held in equivalent traps, and as such we can draw upon 
the relevant literature. In modelling the rotational d 5 mamics for systems of this type, it is a well known problem 
that the equations of motion cannot be solved using spherical polar coordinates iS6ll . Working with the quaternion 
formalism Issll . we built a numerical simulation in XMDS is3l which modelled the coupled dynamics of the rotary 
and COM motion, and experimented with the effects of different dipole moment magnitudes. We expect that packing 
the maximum possible charge onto the surface of the sphere would naturally lead to the most isotropic distribution 
isih . which in our case gives us q = 1500 eV. We find that the effect of the rotary motion does cause an appreciable 
heating effect for a large enough dipole moment, which occurs when qo = R x 5 eV. This limit on the dipole 
magnitude qo is perhaps more clearly expressed as a ratio gedge/fZcore < 1/300, where we have simplified the charge 
distribution into q = qcore + 9edge/ 9core being the charge which is completely isotropic and effectively bunched at 
the COM, and gedge being the charge which is bunched at one position of the surface, and hence causes the dipole 
moment. 

We can infer a bound on the ratio from experimental data already gathered. As stated, this experimental 
proposal rests upon the hybrid type trap which has already been constructed and used. In the experiments which 
we have conducted so far, we have used far lower charge numbers, which would presumably distribute themselves 
more anisotropically the situation presented in this paper. In analysing the dynamics of these lower-charge systems 
we do not find strong evidence of the coupling between the rotational motion to the translational dynamics of the 
sphere. We take this as an indication that in these situations we indeed have Therefore, we conclude that 

in our scenario, in which we would presumably have a more isotropic distribution, we would have an even lower 
ratio and therefore an acceptably low heating rate from the rotational dynamics. 
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